42 research outputs found
Planewave density interpolation methods for 3D Helmholtz boundary integral equations
This paper introduces planewave density interpolation methods for the
regularization of weakly singular, strongly singular, hypersingular and nearly
singular integral kernels present in 3D Helmholtz surface layer potentials and
associated integral operators. Relying on Green's third identity and pointwise
interpolation of density functions in the form of planewaves, these methods
allow layer potentials and integral operators to be expressed in terms of
integrand functions that remain smooth (at least bounded) regardless the
location of the target point relative to the surface sources. Common
challenging integrals that arise in both Nystr\"om and boundary element
discretization of boundary integral equation, can then be numerically evaluated
by standard quadrature rules that are irrespective of the kernel singularity.
Closed-form and purely numerical planewave density interpolation procedures are
presented in this paper, which are used in conjunction with Chebyshev-based
Nystr\"om and Galerkin boundary element methods. A variety of numerical
examples---including problems of acoustic scattering involving multiple
touching and even intersecting obstacles, demonstrate the capabilities of the
proposed technique
Harmonic density interpolation methods for high-order evaluation of Laplace layer potentials in 2D and 3D
We present an effective harmonic density interpolation method for the
numerical evaluation of singular and nearly singular Laplace boundary integral
operators and layer potentials in two and three spatial dimensions. The method
relies on the use of Green's third identity and local Taylor-like
interpolations of density functions in terms of harmonic polynomials. The
proposed technique effectively regularizes the singularities present in
boundary integral operators and layer potentials, and recasts the latter in
terms of integrands that are bounded or even more regular, depending on the
order of the density interpolation. The resulting boundary integrals can then
be easily, accurately, and inexpensively evaluated by means of standard
quadrature rules. A variety of numerical examples demonstrate the effectiveness
of the technique when used in conjunction with the classical trapezoidal rule
(to integrate over smooth curves) in two-dimensions, and with a Chebyshev-type
quadrature rule (to integrate over surfaces given as unions of non-overlapping
quadrilateral patches) in three-dimensions
Windowed Green Function method for layered-media scattering
This paper introduces a new Windowed Green Function (WGF) method for the
numerical integral-equation solution of problems of electromagnetic scattering
by obstacles in presence of dielectric or conducting half-planes. The WGF
method, which is based on use of smooth windowing functions and integral
kernels that can be expressed directly in terms of the free-space Green
function, does not require evaluation of expensive Sommerfeld integrals. The
proposed approach is fast, accurate, flexible and easy to implement. In
particular, straightforward modifications of existing (accelerated or
unaccelerated) solvers suffice to incorporate the WGF capability. The
mathematical basis of the method is simple: the method relies on a certain
integral equation posed on the union of the boundary of the obstacle and a
small flat section of the interface between the penetrable media. Numerical
experiments demonstrate that both the near- and far-field errors resulting from
the proposed approach decrease faster than any negative power of the window
size. In the examples considered in this paper the proposed method is up to
thousands of times faster, for a given accuracy, than a corresponding method
based on the layer-Green-function.Comment: 17 page
Boundary integral equation methods for the solution of scattering and transmission 2D elastodynamic problems
We introduce and analyse various regularized combined field integral equations (CFIER) formulations of time-harmonic Navier equations in media with piece-wise constant material properties. These formulations can be derived systematically starting from suitable coercive approximations of Dirichlet-to-Neumann operators (DtN), and we present a periodic pseudodifferential calculus framework within which the well posedness of CIER formulations can be established. We also use the DtN approximations to derive and analyse OS methods for the solution of elastodynamics transmission problems. The pseudodifferential calculus we develop in this paper relies on careful singularity splittings of the kernels of Navier boundary integral operators, which is also the basis of high-order Nystrom quadratures for their discretizations. Based on these high-order discretizations we investigate the rate of convergence of iterative solvers applied to CFIER and OS formulations of scattering and transmission problems. We present a variety of numerical results that illustrate that the CFIER methodology leads to important computational savings over the classical CFIE one, whenever iterative solvers are used for the solution of the ensuing discretized boundary integral equations. Finally, we show that the OS methods are competitive in the high-frequency high-contrast regime.Catalin Turc gratefully acknowledges support from National Science Foundation (NSF) through contract DMS-1614270 and DMS-1908602